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^ , Abstract: Although imphcit in the discovery of the Schwarzschild solu- 

tion 40 years earlier, the issues raised by the theory of what are now known 



O 



O ■ as black holes were so unsettling to physicists of Einstein's generation that 

^ ' the subject remained in a state of semiclandestine gestation until his demise. 

■*»>> . That turning point - just half a century after Einstein's original foundation 

qh| of relativity theory, and just half a century ago today - can be considered 



^ 



to mark the birth of black hole theory as a subject of systematic develop- 
PfJ' ment by physicists of a new and less inhibited generation, whose enthusastic 

investigations have revealed structures of unforeseen mathematical beauty, 
^ ■ even though questions about the physical significance of the concomitant 

singularities remain controversial. 
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Figure 1: Numerically simulated view of (isolated) spherical black hole illu- 
minated only by uniform distant sky background from which light is received 
only for viewing angle a > (3 where, for an observer falling radially from 
rest at large distance, the angle (5 subtended by the hole will be given by 
the formula (J45|l obtained in the appendix. 



1. Introduction: Schwarzschild's unwelcome solution. 

This illustrated review is intended to provide a brief overview of the emer- 
gence, during the last half century, of the theory of ordinary (macroscopic 
4- dimensional) black holes, considered as a phenomenon that (unlike the 
time reversed phenomenon of white holes) is manifestly of astrophysical im- 
portance in the real world. The scope of this review therefore does not cover 
quantum aspects such as the Bekenstein-Hawking particle creation effect, 
which is far too weak to be significant for the macroscopic black holes that 
are believed to actually exist in the observable universe. Nor does it cover 
the interesting mathematics of higher dimensional generalisations, a subject 
that is (for the time being) so far from relevance to the known physical world 
(in which - according to the second law of thermodynamics - the distinction 
between past and future actually matters) that its practitioners have formed 



a subculture in which the senior members seem to have forgotten (and their 
juniors seem never to have been aware of have been aware of) the distinction 
between black and white holes, as they have adopted a regretably misleading 
terminology whereby the adjective "black" is abusively applied to any brane 
system that is hollow - including the case of an ordinary (black or white) 
hole, which, to be systematic, should be classified as a (black or white) hollow 
zero brane of co dimension 3. 

The rapid general acceptance of the reality and importance of the positrons 
whose existence was implied by Dirac's 1928 theory of the electron is in strik- 
ing contrast with the widespread resistance to recognition of the reality and 
importance of the black holes whose existence was implied by Einstein's 1915 
theory of gravity. It is symptomatic that black holes were not even named as 
such until more than half a century later. The sloth with which the subject 
has been developed over the years is illustrated by the fact that although 
the simplest black hole solution was already discovered (by Schwarzschild) 
in 1916, the simulation in Figure ^ of the present review (80 years later) 
provides what seems to the first serious reply to the very easy question of 
what it would actually look like, all by itself, with no illumination other than 
that from a uniform sky background. 

Much of the responsibility for the delay in the investigation of the conse- 
quences of his own theory is attributable P] to Einstein himself. Although his 
work had revolutionary implications, Einstein's instincts tended to be rather 
conservative. It was as a matter of necessity (to provide an adequate account 
first of electromagnetism and then of gravitation) rather than preference that 
Einstein introduced the radically new paradigms involved first in his theory 
of special relativity, just a hundred years ago, and then in the work on general 
relativity that came to fruition ten years later. When cherished prejudices 
were undermined by the consequences, Einstein was as much upset as any of 
his contempory colleagues. It could have been said of Albert Einstein (as it 
was said of his illustrious and like minded contempory, Arthur Eddington) 
that he was always profound, but sometimes profoundly wrong. 

The most flagrant example was occasionned by Friedmann's prescient 
1922 discovery of what is now known as the "big bang" solution of the general 
relativity equations, which Einstein refused to accept because it conflicted 
with his unreasonable prejudice in favor of a cosmological scenario that would 
be not only homogeneous (as actually suggested by subsequently available 
data) but also static (as commonly supposed by earlier generations) despite 
the incompatibility (in thermal disequilibrium) of these alternative simpli- 




Figure 2: John Wheeler with Robert Dicke (Princeton 1971). 

fications with each other and with the obvious observational consideration 
(known in cosmologically minded circles as the Cheseaux-Olbers paradox) 
that - between the stars - the night sky is dark. Einstein's incoherent atti- 
tude (reminiscent of the murder suspect who claimed to have an alibi as well 
as the excuse of having acted in self defense) lead him not only to tamper with 
his own gravitation equations by inclusion of the cosmo logical constant, but 
anyway to presume without checking that Friedmann's (actually quite valid) 
solution of the original version must have been mathematically erroneous. 

Compared with his tendency to obstruct progress in cosmology, Einstein's 
conservatism was rather more excusible in the not so simple case of what are 
now known as black holes. It is understandable that (like Eddington) he 
should have been unwilling to explore the limitations on the validity of his 
theory that are indicated by the weird and singular - or as Thorne pQ puts it 
"outrageous" - features that emerge when strong field solutions of the general 
relativity equations are extrapolated too far into the non linear regime. 

At the outset Einstein's interest in the spherical vacuum solution of his 
1915 gravitational field equations was entirely restricted to the weak field 



regime, far outside the "horizon" at ?' = 2m in the simple exact solution 

ds^ = r2(d^2 ^ sin2^^^2^) ^ drV(l-2m/r) - (l-2m/r)dt^ (1) 

that was obtained within a year, but that was immediately orphanned by 
the premature death of its discoverer, Karl Schwarszschild, after which its 
embarrassing physical implications were hardly taken seriously by anyone - 
with the notable exception of Oppenheimer |2] - until the topic was taken 
up by a less inhibited generation subsequent to the death of Einstein himself, 
just half a century ago, at Princeton in 1955. It was only then (and there) 
that John Wheeler inaugurated the systematic development of the subject ~ 
for which he coined the name "black hole" theory - in a series of pionnering 
investigations that started |3] by addressing the crucial question of stability, 
while not long afterwards, on the other side of the "iron curtain" another 
nuclear arms veteran, Yacob Zel'dovich, initiated an independent approach 
in to the same problem (using the alternative name "frozen star" which in 
the end did not catch on). 

2. Outcome of stellar evolution: Chandra's unwelcome limit. 

The question of gravitational trapping of light had been raised in the 
eighteenth century by Michel and Laplace, whose critical mass m ^ p~^/^ 
assumed the standard mass density that is understood (on the basis of quan- 
tum theory as developped by 1930) to result in hadronic matter from balance 
between Fermi repulsion and electrostatic attraction which (in Planck units, 
with proton and electron masses nip ^ 10^^^ , m,e ~ 10~^^ gives p ~ nipn 
with n f» A~^ for the Bohr radius A ~ e^/me with e^ ~ 1/137 . However 
most theorists refused to face the issue of gravitational collapse even after 
progress in quantum theory lead to Chandrasekhar's 1931 discovery of the 
maximum mass m ?a m" ^ for cold body - which is attained when relativis- 
tic gas pressure P ~ n/X ^ n'^^^ provides the support required by virial 
condition P ^ nip p^l'^ . 

For a lower mass m < m~'^ , stellar evolution at finite temperature 0, 
with gas pressure P ~ n© subject to p ~ Q^/mprn^ , can terminate in 
cold equilibrium supported by non-relativistic Fermi pressure P ~ n^^^ /rrie 
giving p Ri m^rrLpm^ for a white dwarf, or P ?a iv'/^ /rrip giving p Ri rripm^ 
for neutron star, as shown in Figure 01 

However a self gravitating mass of hot gas will be radiation dominated 
with P ^ Q^ , whenever its mass exceeds the Chandrasekhar limit, m > 
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Figure 3: Logarithmic plot of density versus mass 
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Figure 4: Conformal Projection diagrams showing firstly the combined black 
and white hole geometry obtained (belatedly jB] in 1960, by Wheeler in col- 
laboration with Kruskal) as the artificial analytic (vacuum) extension of the 
Schwarzschild solution, and secondly a more astrophysically natural exten- 
sion with homogeneous interior (found in 1939 [2j by Oppenheimer in col- 
laboration with Sneyder) in which the shaded egion is the "domain of outer 
communications" and the unshaded region is the prototype example of a hole 
qualifiable as black in the strict sense. 



m~ , so that, as first understood by Chandra's Cambridge research director, 
Arthur Eddington, its condition for (thermally supported) equilibrium will be 
given by p ~ 0^/m-^/^. What Chandra could never get Eddington to accept is 
that, for such a large mass, no cold equilibrium state will be available, so after 
exhaustion of fuel for thermonuclear burning (at 6 ~ e^mp ) gravitational 
collapse will become inevitable. 

3. Spherical collapse past the horizon 

Eddington's example shows how, as has described in detail by Werner 
Israel [S] (and in striking contrast with the open mindedness of Michel and 
Laplace a century and a half earlier) physicists of Einstein's generation tried 
to convince themselves that nature would never allow compacification within 
a radius comparable to the Schwarzschild value. While Einstein lived, even 
after Chandrasekhar's discovery had shown that such a fate might often be 
difficult to avoid, the implications were taken seriously only by Oppenheimer 
and his colleagues, who showed |2] how, as shown in Figure |3 the solutions 
of Schwarzschild and Friedmann could be combined to provide a complete 
description of the collapse of a homogeneous spherical body through what is 
now called its event horizon all the way to a terminal singularity. 

Despite the persuasion of such experienced physicists as Wheeler and 




Figure 5: First realistic simulation of distant view of spherical black hole 
with thin accretion disc by Jean-Pierre Luminet, 1978. 

Zel'dovich, and the mathematical progress due to younger geometers such as 
Robert Boyer and particularly Roger Penrose, the astrophysical relevance of 
the region near and within the horizon continued to be widely disbelieved 
until (and even after) the 1967 discovery [7j by Israel of the uniqueness of 
the Schwarzschild geometry as a static solution: many people (for a while 
including Israel himself jS]) still supposed (wrongly) that the horizon was an 
unstable artefact of exact spherical symmetry. It is therefore not surprising 
that the question of what such a black hole would actually look like was not 
addressed until much more recently, particularly considering that nothing 
would be seen at all without some source of illumination. 

The realisation that many spectacular astrophysical phenomena ranging 
in scale from supermassive quasars in distant parts of the universe down 
to stellar mass X ray sources within our own galaxy may be attributed to 
accretion discs jHl El E] round more or less massive black holes has how- 
ever provided the motivation for increasingly realistic numerical simulations 
(Figures El and inj of what would be seen from outside in the presence of an 
illuminating source of this kind ^3 El • 

As the most easily calculable example, I have shown in the appendix how 
to work out the case shown in Figure ^ of an isolated spherical black hole for 
which the only source of illumination is a uniform distant sky background, 




Figure 6: Simulation of close up view of spherical black hole with thin accre- 
tion disc by Jean Alain Marck, 1996. 



viewed as a function of proper time, 



Am / r \3/2 



(2) 



by a (doomed) observer falling towards the singularity inside the black hole, 
with zero energy and angular momentum. In such a case the redshift Z 
determining the observed energy S/{1 + Z) of a photon emitted from the 
sky background with the uniform average energy S say will be given by the 

formula 

cos a 
Z = ^===, a>/5, (3) 

where a is the apparent angle of reception, which must of course excede 
the apparent angle fi subtended by the black hole. This means that the 
redshift will be positive (so that the sky will appear darker than normal) 
due to the Doppler effect, for photons coming in from behind the observer 
(with a > Tc/2 ). However photons received in the range j3 < a < 7r/2 will 
be blueshifted by an amount that will diverge, as shown in Figure [71 as the 
singularity is approached. 
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Figure 7: Plot of reception angle a against countdown proper time r for 
arrival at singularity of radially falling zero-energy observer in units such 
that 2'ni = ■\/2/3 so that null orbit radius r = 3m is crossed when r = — 1. 
Constant brightness contours indicate intensity of light received from uniform 
background, which will be inversely proportional to 4th power of redshift 
factor {1 + Z) , having unshifted value for rays arriving at angle n/2 . 



4. Discovery of horizon stability and of Kerr solution 

Following the demise of Einstein (and the development of nuclear weapons) 
a new (less inhibited) generation of physicists, lead by Wheeler and Zel'dovich, 
came to recognise the likelihood - and need in any case for testing - of stabil- 
ity with respect to non-spherical perturbations of what was termed a "black 
hole". Work by Vishweshwara ^2]; Price ^1], and others confirmed that 
"anything that can be radiated away will be radiated away" - leaving a fi- 
nal equilibrium state characterised only by mass and angular momentum. 
The (still open) mathematical question of the extent to which this remains 
true (with singularities hidden inside horizon) for very large deviations from 
sphericity was raised by the "cosmic censorship" conjecture formulated by 
Roger Penrose ^Hl ^j but in any case the relevance of black holes for astro- 
physical phenomena (notably quasars) was generally accepted in astronomi- 
cal circles from 1970 onwards. 

The generic form of what was afterwards recognised to be the final black 
hole equilibrium state state in question was discovered in 1963, when Roy 
Kerr announced f?! E] that "among the solutions . . . there is one which is 
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Figure 8: Sketch showing projections of null cones in equatorial space section 
of Kerr metric. The pale curve marks the "ergosurface" bounding the region 
where stationary motion is allowed, and the heavy black curve marks the 
horizon bounding the trapped "black hole" region. 

stationary ... and also axisymmetric. Like the Schwarzschild metric, which 
it contains, it is type D ... m is a real constant ... The metric is 

ds^ = (r^ + a^cos'^d){dd^ + sin^^d^^) + 2(dn + asin^O d(j)){dr + asm'^dd(j)) 

1 5 5 ^ (dM + asin2^d0)^ (4) 

r^^ + a^cos^py 

where a is a real constant. This may be transformed to an asymptotically 
flat coordinate system ... we find that m is the Schwarzschild mass and ma 

the angular momentum " . 

Since the black hole concept had still not been clearly formulated then, 
it was at first (wrongly) supposed that the physical relevance of this vacuum 
solution would be as the exterior to a compact self gravitating body like a 
neutron star, as suggested by Kerr's (off the mark) conclusion jTTj that it 
would be "desirable to calculate an interior solution." 

What actually makes the Kerr metric so important however, as can be 
see from Figure IHl (using CP. diagrams, which were originally developed for 
this purpose) is the feature first clearly recognised fOl 1201 by Bob Boyer in 
1965, which is that for a^ < m^ the distant sky limit known as "asymptopia" 
is both visible and accessible only in a non-singular "domain of outer com- 
munications" bounded by past and future null (outer) horizons, on which 
A = 0, where 

r = m-|-c, c = yrn? — a? (5) 
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The topology within the black (and white) hole regions was first eluci- 
dated J2I] in terms of Conformal Projections on the symmetry axis in 1966 
and then completely |^ |2Sl 121] by Boyer, Lindquist and myself in 1967 
and 1968 - the year when the much needed term "black hole" was finally 
introduced by Wheeler to describe the region from which light cannot escape 
to "asymptopia" . (A "white hole" region would be one that could not re- 
ceive light from "asymptopia".) In the generic rotating case (unlike the static 
Schwarzschild limit) the well behaved domain outside the black hole horizon 
includes an "ergosphere" region where, as shown in Figure |Hl the Killing vec- 
tor generating the stationarity symmetry becomes spacelike, so that (globally 
defined) particle energies can be negative. 

In contrast with the good behavior of the outer region, r > m + c, I found 
that, as well as having the irremovable ring shaped curvature singularity 
already noticed by Kerr where r^ -|- a^cos^6' — »• 0, , the inner parts of the 
rotating Kerr solutions would always be causally pathological, due to the 
existence near the ring singularity of a small region (see Figure irUjl where 
the axial symmetry generating Killing vector becomes timelike PHll25j . This 
feature gives rise to a causality violating "time machine region" (a feature so 
"outrageous" as to be unmentionable even by Thorne |Tj) that would extend 
all the way out to "asymptopia" (meaning r -^ +oo ) in the - presumably 
unphysical - case for which a^ > ni^ . (I would emphasize that this kind of 
time machine, like those recently considered by Ori [22], would survive even 
if one takes the covering space, unlike a time machine of the wormhole kind 
discussed by Thorne P which is merely an artefact of multiply connected 
space time topology). 

In so far as the (physically relevant) black hole cases characterised by 
a^ < irv^ are concerned, the good news [SHI 12^] (for believers in causality) 
is that the closed timelike lines are all contained within the inner region 
r < m — c. The boundary of the time machine region is constituted by the 
"inner horizon" , where r = m — c. which acts as a Cauchy hypersuface from 
the point of view of inital data for formation of the black hole by gravitational 
collapse. Unlike the outer horizon r = ni + c, whose stability throughout 
the allowed range < a^ < m^ has been even confirmed by Whiting [211, i^ 
was to be expected J2Z1 I2H1 that a Cauchy horizon of the kind occurring at 
r = m — c would be unstable, and it has been shown that outcome is likely 
to be the formation of a curvature singularity of the weak kind designated 
by the term "mass inflation" pTH 11101 m] . 
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5. Seductive mathematical features of Kerr type metrics 

In his original 1963 letter ^7], and with Alfred Schild ^H] in a sequel, 
Kerr obtained the useful alternative form 

ds^ = g^y dx^dx'" = rj^y + 2{m/U)n^ny (6) 

with null covector n^dx^ = du + asin^^ d0 , for f/ = (r^ + a^cos^^)/?' , in 
a flat background. The latter was obtained in the Minkowski form, 

7]^y dx^dx^ = dx^ + df + dz^ - dt^ , (7) 

by setting t = u — r , z = a cos 9 , x + iy = {r — m)e"^sin 9 , which gave 

, „ ,- zdz (rx — ay)dx+(ry+ax)dy , , 

n^dx'^=dt+ +^ 2,2 -^- 8 

r r^ + a^ 

(This form of pure vacuum solution was generalised to higher dimensions by 
Myers and Perry [32] • It is perhaps of greater current cosmo logical interest 

- in view of the evidence that the expansion of the universe is accelerating 

- that this form has also beeen extended to include a cosmological constant 
in a 4 dimensional De Sitter background by myself [23 Ell, while further 
generalisations to a De Sitter background in 5 and higher dimensions p^ IHH] 
have been obtained more recently.) 

As well as time and axial symmmetry, the Kerr solution has a discrete 
PT symmetry that was predictable from Papapetrou's "circularity" theorem 
[HZI, and made manifest in 1967 [22] by the Boyer Lindquist transformation 

dt = dn - (r^ + a^) A"Mr , dip = -dcp + a A"Mr , (9) 

with 

A = r^ - 2mr + a^ . (10) 

This gives Kerr's null form as 



n^dx'^ = dt - a sin^^ d(/7 + ^^ A M?' , p= Vr^ + o^cos^^, . (11) 

The metric itself is thereby obtained in the convenient form 

ds^ = g^(^ + d9A+{r^W)sm^9dip^ + '^{dt-asin'9dipy-dt^ , (12) 
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Figure 9: Representation of equatorial space-time section of Kerr metric by 
Conformal Projection diagrams. (Such C.P. diagrams were originally devel- 
oped for this purpose.) The second version achieves complete compactifiction 
by letting scale for successive universes tend to zero at extremities of chain. 



in which there are cross terms involving the non-ignorable differentials, dr 
and d6 , but - as the price for this simplification - if a^ < m^ there will be 
a removable coordinate singularity on the null "horizon" where A vanishes. 
Whereas the possibility of making the foregoing simplification was pre- 
dictable in advance, there was no reason to anticipate the discovery [221 EH] 
that, in addition to the ordinary "circular" symmetry generated by Killing 
vectors, k^d/dx'^ = d/dt and h^d/dx^ = d/dip , the Kerr metric would 
turn out to have the hidden symmetry that is embodied in the canonical 
tetrad 



9,^=J2^A-K^ 



(13) 



specified by 



/dx'^ = (^/v^)dr, 



rdx^: 



gdO, 



(14) 



'dldxf' {r'^ + a'^)dLp-adt ^"dx'^ dt-asin^Odip 



sm( 



A 



(15) 
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In terms of this canonical tetrad, the Kerr-Schild form of the metric is 
expressible as 

9,u = r/^. + 2mr(<+<)(<+<) , (16) 

while the Killing- Yano 2-form brought to light by Roger Penrose and his 
coworkers is expressible as 

f,, = 2acose^\/^^-2r^\Xy (17) 

The property of being a Killing- Yano 2-form means that it is such as to 
satisfy the very restrictive condition condition 

V^/^p = V[^/vp] , (18) 

thus providing a symmetric solution 

K,. = f,,n, V(^ir,p) = 0, (19) 

of the Eisenhart type Killing tensor equation, as well as secondary and 
primary solutions kf" = K^'^k'' =a^k>' + ah>' and k^" = ke^^^^Vyfpa of the 
ordinary Killing vector equation V(^A;j,) = . 

For affine geodesic motion, jfV^p^ = 0, one thus obtains (energy and axial 
angular momentum) constants S = k^p^ and A4 = h^Pi,, while the Killing 
tensor gives the constant /C = K^^^p^p^ = i^i^ , with (angular momentum) 
J^ = f,.uP^ obeying p'^VJ^ = 0. 

There will also [lOj be corresponding (self adjoint) operators 

S = ik''Vu, M=ih''Vu, IC = Vf.K'^'V, , (20) 

whose action on a scalar field commute with that of the the Dalembertian 
D = VVu '■ in other words [£,0] = 0, [A^, D] = 0, and (consistently with 
the integrability condition KP\^fj_Ry^p = ) also [/C, D] = . 

The ensuing integrability of the geodesic equation |2S1 and of the scalar 
wave equation is equivalent to their solubility by separation of variables |23| 
EH] • The possibility of extending these rather miraculous separability prop- 
erties to the neutrino equation jll] and even to the massive spin 1/2 field 
|12l Eni as governed by the Dirac operator V = 7'^V^ is attributable to 
corresponding spinor operator conservation laws 

[S,V] = 0, [M,V] = 0, [J,V] = 0, (21) 
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Figure 10: Reproduction from 1972 Les Houches notes j21]: sketch of {r,9} 
section through ring singularity at junction between heavily shaded region 
responsible for causality violation where axisymmetry generator is spacelike, 
and lightly shaded ergo region where time summetry generator is spacelike. 

of energy, axial angular momentum, and (unsquared) total angular momen- 
tum, as respectively given jUj by 



S = ik''V. + \z{V,K)rY , M=ih'V, + \z{V^K)YY , (22) 

and 

j = ^nYf,'''^u-K). (23) 

Such a neat commutation formulation is not (yet?) available for Teukol- 
sky's extension [13 HEj of solubility by separation of variables to massless 
spin 1 and spin 2 fields representing electromagnetic and gravitation pertur- 
bations - of which the latter are particularly important for Bernard Whiting's 
demonstration j^H] of stability. An even more difficult problem is posed by 
the charged generalisation |32| of the Kerr black hole metric, which retains 
many of its convenient properties (and is noteworthy for having the same 
gyromagnetic ratio as the Dirac electron [221 EHl) but which gives rise to a 
system of coupled electromagnetic and gravitational perturbations that has 
so far been found to be entirely intractible. 
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6. No hair and uniqueness theorems for black hole equilibrium 

The overwhelming importance of Kerr solution derives from its provision 
of the generic representation of the final outcome of gravitational collapse, 
as was made fairly clear in 1971 by the prototype no-hair theorem [^IHU] 
proving that no other vacuum black hole equilibrium state can be obtained by 
continuous axisymmetric variation from the spherical Schwarzschild solution 
that had been shown by the earlier work of Israel [7j (before the generic 
definition of a black hole was available) to be only static possibility. 

Conceivable loopholes (such as doubts about the axisymmetry assump- 
tion) in the reasonning leading to this conclusion (which was rapidly - per- 
haps too uncritically - accepted in astronomical circles) were successively 
dealt with by the subsequent mathematical work of Stephen Hawking J3I], 
David Robinson [S2] and other more recent contributors |^ |33J EH] to what 
has by now become a rather complete and watertight uniqueness theorem 
for pure vaccum black hole solutions in 4 spacetime dimensions. It should 
however be remarked jSH] that there are some mathematical loose ends (con- 
cerning assumptions of analyticity and causality) that still need to be tidied 
up. 

The demonstration uses ellipsoidal coordinates for the 2-dimensonal space 
metric ds^ = dA^/(A^ — c^) -t- d/.i^/(l — /.i^), in terms of which the generic 
stationary axisymmetric asymptotically fiat vacuum metric is known from 
the work of Papapetrou ^7\ to be expressible in the form 

ds^ = g^ds^ + Xidif -udtf - {\^-c^){l-iJ^)dt^ . (24) 

for which, by the introduction of an Ernst [SZj type potential given by 
X'^duj/dX = {l—fi'^)dY/dfi, the relevant Einstein equations will be obtainable 
from the (positive definite) action /dAd/i(|VXp + \VY\'^)/X'^. 

The black hole equilibrium problem is thus [IHl 1211 HO] reduced to a non 
linear 2 dimensional elliptic boundary value problem for the scalars X, Y, 
subject to conditions of regularity on the horizon (with rigid angular velocity 
Q ) where A = c and to appropriate boundary conditions on the axis where 
fi = ±1 and at large radius A — >■ cxd in terms of angular momentum ma . 

The uniqueness theorem states that this 2 dimensional boundary problem 
has no solutions other that those given given (with X = r — m, fi = cos 9) by 
the Kerr solution having mass ??7 = y/c^+a? and horizon angular velocity 
Q = a/2m{m+c) . The proof is obtained from an identity equating a quantity 
that is a positive definite function of the relevant deviation (of some other 
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caustic 



Figure 11: Sketch (in plot of time against space) of segment of black hole 
horizon showing how a null generator (obtained as limit of escaping timelike 
curves) can begin (on a caustic) but can never end towards the future. 



hypothetical solution from the Kerr value) to a divergence whose surface 
integral can be seen to vanish by the boundary conditions. 

The original no-hair theorem (applying just to the small deviation limit) 
was based on an infinitesimal divergence identity that I obtained by a hit and 
miss method [IH] that was generalised by Robinson [S2] to the finite differ- 
ence divergence identity that was needed to complete the proof in the pure 
vacuum case. For the electromagnetic (Einstein Maxwell) generalisation, the 
analogous step from an infinitesimal no-hair theorem |58j to a fully non-linear 
uniquenes theorem was more difficult, and was not obtained until our hit and 
miss approach was superceded by the more sophisticated methods that were 
developed later on by Mazur jSHl ED] and Bunting jHH I 



6. Further developments 

After it had become clear that (in the framework of Einstein's theory) 
the Kerr solutions (with a < m ) are the only vacuumm black hole equilib- 
rium states, the next thing to be investigated was the way the black holes 
will evolve when the equilibrium is perturbed. A particularly noteworthy 
result, based on concepts (see Figure [TT] ) developed in collaboration with 
Penrose jHl] was the demonstration by Stephen Hawking jMJ EI] that the 



area of a black hole horizon (which is proportional to what Christodoulou 
^^ had previously identified as irreducible mass) can never decrease. More 
particularly it was shown jBH] that the area would grow, not only when the 
hole swallowed matter but more generally whenever the null generators of 
the horizon were subjected to shear. It was remarked that this effect could 
be described in terms of an effective viscosity and that the horizon could also 
be characterised jHIi IHH] by an effective resistivity. 

Later astrophysical developments were concerned more with surrounding 
or infalling matter - for example in accretion discs - than with the black 
hole as such, at least until recently. However the prospect of detecting grav- 
itational radiation in the foreseeable future has encouraged a resurgence of 
interest in purely gravitational effects, particularly those involved in binary 
coalescence. The climax of a coalescence is too complicated to be dealt with 
except by advanced methods of numerical computation, but the quasi sta- 
tionary preliminary stages are more amenable jHSl HOI IHl IZ21 ESI Ell , ^^ also 
are the final stages of ringdown, which can be analysed in terms of quasi 
normal modes (and their superpositions in power law tails of the kind first 
described by Price ^3]) which have been the subject of considerable atten- 
tion, particularly concerning the infiuence of rotation [75117^17711751 17^. 

Appendix: null geodesies in spherical case. 

Although it would be insufficient for the complete Kruskal (black and 
white hole) extension, in order to cover a purely black hole (Oppenheimer 
Sneyder type) extension of the Schwarzschild solution, it will suffice to use 
an outgoing null coordinate patch of the kind introduced for the Kerr metric 
(HI) for which, when a = 0, the metric will be given in terms of x" = u, 
x^ = r, x"^ =9, x^ = (p, , simply by 

ds^ = -{1 -2m/r) du^ + 2dudr + r'^de'^ + sin^e d^^ . (25) 

Within such a system, an observer falling in freely from a large distance 
with zero energy and angular momentum will have a geodesic trajectory 
characterised by fixed values of the angle coordinates 6 and (p and by a 
radial coordinate r that is given implicitly as a monotonically decreasing 
function of proper time by (0) and as a monotonically decreasing function of 
the ignorable coordinate u by the relation 



^^^^^ = -a/— f— + 3)- — -21n|l + J^^ (26) 

2m 3\2m\2m J 2m \ \ 2m ' ^ ^ 
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(0) 



in which uq is a constant of integration specifying the value of u for which 
the trajectory terminates at the singular limit r — >• . For such a trajectory 
the (future oriented) timelike unit tangent vector will be given by 

[l + ,/2^y\ e^^^ = -y/2^, (27) 

and the tetrad specifying a corresponding local reference frame can be com- 
pleted in a natural manner by using the associated (outward oriented) or- 
thogonal spacelike unit vector, which will be given by 

e^» =(l + v/2W^)"\ e^^^^ = l, (28) 

together with two other horizontally oriented unit vectors whose specification 
will not matter for the present purpose because of the rotation symmetry of 
the system. 

Let us consider the observation of photon that arrives with trajectory 
deviating by an angle a say from the outward radial direction. The first two 
components of its (null) momentum vector, p^^ say, will evidently be given 
in terms of its energy S with respect to such a frame by 

j/ = ^(e;^+cos«e;j, (29) 

for i = 0,1 . This locally observed energy is to be compared with the 
globally defined photon energy (as calibrated with respect to the asympotic 
rest frame at large distance) that will be given in terms of the timelike Killing 
vector with components k^ = 5^ by 

S = -k^if = (1 -2m/r) p" -p^ (30) 

the important feature of the latter being that it is conserved by the affine 
transport of the momentum vector along the null geodesic photon trajectory 
to which it is tangent. It can thus be seen that the corresponding locally 
observed energy S will be related to the globally defined energy constant S 

by 

£ = £il + Z), (31) 

with the redshift Z given by Q, and that the associated component ratio 
will be given by 



P ^ , I T^ I cos a — \/2mh \ 

— = 1 + J2m r) y- '— . 32 

po ^ V / ; I cosa + 1 ' ^ ' 



20 



As the (unsurprising) spherical hmit of the (still rather mysterious) sepa- 
rability of Kerr's rotating generalisation, the evolution of the relevant affinely 
transported momentum components will be given in terms of the energy con- 
stant £ and the associated squared angular momentum constant /C [22] 
(using a dot for differentiation with respect to the affine parameter) by 

„ . ±Vtz + v ^ . ±Vn ,„„. 

p° =u = — — , pi = r = — — , (33) 

r^ — Zmr r^ 

where 

V = Er^ , 7^ = p2 - /C(r2 - 2mr) . (34) 



We thereby obtain 



£!^|^^,:-W.r.(:.-^). ,35) 



and hence, by comparison with 



\fll, /— / cos a — \J'lmlr \ 

V \ \/r j2m — cos a j 

This expression can be used to evaluate the squared angular momentum 
constant as a function of the locally defined energy 8 and angle a in the 
form 

K = r^^'sin^a (37) 

in which the variable £ will itself be given via ()31|) in terms of the globally 
defined energy constant £ by the - red or blue - shift formula Q so that for 
the square of the constant ratio of angular momentum to energy one obtains 

^ ^ ( rsina \ 

£^ \l-^2^-cosa) ■ ^ ' 

With respect to an unconventional affine parameter orientation condition 
to the effect that the energy should always be non-negative, iS" > , it can be 
seen (in view of the consideration that the squared angular momentum con- 
stant must necessarily be non-negative, /C > ) that a null geodesic segment 
will be appropriately be describable as "incoming" or "outgoing" according 
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to whether the upper or lower of the sign possibihties ± is apphcable, i.e. ac- 
cording to whether the right hand side of (|36|) is positive or strictly negative. 
It is however to be remarked that this convention will be consistent with the 
usual requirement that the affine parameter orientation be future oriented, 
giving u > , only outside the horizon and for "ingoing" null segments within 
the horizon, where r < 2m , but that for "outgoing" null segments within 
the horizon it would entail the opposite orientation convention, giving n < 
. With respect to the usual parameter orientation condition giving it > the 
"outgoing" null segments within the horizon will need to be parametrised the 
other way round, which means that they will be characterised by negative 
energy S < and by the upper of the sign possibilities ± . 

Whichever convention is used, it can be seen that within the horizon the 
radius r will aways be a decreasing function of u , even for the (relatively) 
"outgoing" null segments, and that only an "incoming" null segment can 
cross the horizon at a finite value of n . It can be seen from (jHUj) that 
outside the horizon (i.e. for r > 2m ) the criterion for a null segment to be 
classified as "incoming" is that it should have cos a < y2rn/r and that the 
corresponding requirement within the horizon is cos a < yr/2m . 

It can be deduced from the expression ()H4j) that the function 7^ will 
remain positive wherever r is positive if IC/S"^ < 27?n^ . In such a case, the 
null geodesic will either be permanently "incoming", proceding all the way 
from "infinity" (i.e. the limit r ^ oo ) down to the internal singularity (i.e. 
the limit ?■ — > ), or else it will be permanently "outgoing", proceding all 
the way to the singularity or to infinity depending on whether it inside or 
outside the finite horizon radius value r = 2m to which it extends in the 
infinite past, i.e. as u —>■ — oo . As well the such "ingoing" and permanently 
"outgoing" possibilities, the critical case 

JC/S^ = 27m^ (39) 

includes also the exceptional possibility of a marginally outgoing - effectively 
"trapped" - null trajectory with fixed radius r = 3m . 

When the angular momentum exceeds this critical value, i.e. if IC/S"^ > 
2'lmP' , there will be a forbidden range r^ < r < ?'+ of values of ?' for which 
7^ < . It can be seen from E^ that relevant limits are explicitly obtainable, 
as the non negative solutions of the cubic equation 

£h-l - }C{r± - 2m) = , (40) 
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in the form 

r± = 2V'^73^cos{^±/3}, (41) 

with 

^^ = nT arcsinv^l - 21S'^rrfi/K: , (42) 

which evidently entails the conditions 2m < r_ < 3m < r+ . 

This means that for a value of IC/S"^ above the critical bound (jH^ the 
possible null trajectories will be classifable as "free" or "trapped" . The "free" 
geodesies are initially "incoming" from "infinity" but become "outgoing" 
after reaching the inner bound at r = r+ so as to remain in the range 
r > r_|_ . The "trapped" geodesies are either permanently "ingoing" within 
the horizon or else are initially "outgoing" from just outside the horizion but 
become "ingoing" after reaching the outer bound r = r- so as to remain 
within the range r < r^ . 

For a position in the range r < 3m the only kinds of "bright" geodesic, 
meaning those coming from the distant sky at "infinity", are of the per- 
manently "incoming" kind characterised by IC/S'^ < 27m^ . whereas for a 
position in the range r < 3m there will also be "bright" geodesies of the 
"free" kind characterised by IC/S"^ > 27m^ . Apart from the special case of 
the circular null geodesies at r = 3m , all the other kinds of null geodesic 
can be classified as "dark" since they can be seen to have emerged from near 
the horizon limit radius r — >■ 2m, in the distant past (the limit u — >■ — oo) 
and so can be interpreted as trajectories of very highly redshifted radiation 
from the infalling matter that be presumed to originally formed the black 
hole whose static final state is under consideration here. 

It can be seen that the ratio IC/S'^ specified as a function of coso; 
by ()38j) will be monotonically increasing in the "incoming" range, i.e. for 
— 1 < cos a < y2m/r where r > 2m and for —1 < cos a- < \/r j2]n where 
r < 2m . At the upper end of this "incoming" range the ratio YijE'^ the 
tends to a maximum that will be finite - with value r"^ j{r — 2??i) - outside 
the horizon, but that will be infinite inside the black hole. The ratio /C/^^ 
will then decrease monotonically for the higher "outgoing" part of the range 
of cos a . 

The critical value (jH^ will be attained for two values of cos a , of which 
the lower one, cos a = X_ say, will be in the "incoming" range, and the 
higher one, cos a = X ^ will be in the "outgoing" range. It can be seen from 
(|H^ that these values will be obtainable as the upper and lower roots of 

^'''\-Xl) = 21m^{\-^2^^X^f , (43) 
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Figure 12: As in Figure [7| but with higher time resolution, showing just 
final stages after observer passes horizon at r = 2m (which occurs when 
T = —(2/3)^/^ ~ 0.54 in units used here) as the black hole's apparent 
angular radius jS increases towards its upper limit 7r/2 when ?' — >■ . 



which will be real and distinct except at r = 3m where they will coincide. The 
range of angles characterising the "bright" geodesies will therefore be given 

by 



-1 < costt < COS/9 



(44) 



(so that /3 will be interpretable as the apparent angular radius of the black 
hole) with a bounding value cos/3 that will be given by cos/:^ = X_ , within 
the radius of the circular null trajectory, i.e. for r < 3m, , while in the outer 
regions for which r > 3m it will be given by cos /3 = X^ . 

The required solutions of ()43|1 are expressible in terms of the dimensionless 
variable f = r/2m by X^ = {27^/f±\2f^ - 3f\Vf^ + 3f) /{Af^ + 27). It 
can thus be seen that (for the freely falling observer) the apparent angular 
size l3 of the black hole - as shown in the simulation of Figure d and as 
plotted against the proper time ^ in Figure [7| and Figure [T21- will be given 
as a function of the dimensionless radial variable f = r/2m by the analytic 

formula 

27v^ + (2r 



cos/3 



If^ 



3f)Vf'^ + 3f 



4?-3 + 27 



(45) 
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